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A quantum positioning system (QPS) is proposed that can provide a user with aU four 
of his space-time coordinates. The user must carry a corner cube reflector, a good clock, 
^^ ■ and have a two-way classical channel of communication with the origin of the reference 

^^ I frame. Four pairs of entangled photons (biphotons) are sent through four interferometers: 

three interferometers are used to determine the user's spatial position, and an additional 
interferometer is used to synchronize the user's clock to coordinate time in the reference 
frame. The spatial positioning part of the QPS is similar to a classical time-of- arrival (TOA) 
system, however, a classical TOA system (such as GPS) must have synchronized clocks that 
keep coordinate time and therefore the clocks must have long-term stability, whereas in the 
QPS only a photon coincidence counter is needed and the clocks need only have short-term 
stability. Several scenarios are considered for a QPS: one is a terrestrial system and another 
is a space-based-system composed of low-Earth orbit (LEO) satellites. Calculations indicate 
^ . that for a space-based system, neglecting atmospheric effects, a position accuracy below the 

1 cm-level is possible for much of the region near the Earth. The QPS may be used as a 
primary system to define a global 4-dimcnsional reference frame. 

(3 PACS numbers: 06.30.Ft, 06., 95.55.Sh, 42.50.Dv 
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INTRODUCTION 



C^ , During the past several years, the Global Positioning System (GPS) has practically become 

a household word. The GPS is a U.S. Department of Defense satellite system that is used by 
both the military and civilians for navigation and time dissemination [ij, |2, |^. Automobile, ship, 
aircraft, and spacecraft use the GPS for navigation. Telephone and computer network systems 
/\ ' that require precise time use the GPS for time synchronization. The GPS is a complex system 

j^ ■ consisting of approximately 24 satellites orbiting the Earth in circular orbits at approximately 4.25 

Earth radii. The GPS is designed so that signals travel line-of-site from satellite to user, and from 
any place on Earth at least four satellites are in view. If a user receives four GPS satellite signals 
simultaneously from four satellites, s = 1,2,3,4, and the satellites' space-time coordinates [ts^rg) 
at time of transmission are known, the user can solve for his unknown space-time coordinates, 
{to-iVo), by solving the four equations [J, |5| 



^ |2 „2(^ + ^2 



c\to-t,f = 0, s = l,..A (1) 



In Eq.(^ we assume that the signals travel on four light cones that are centered at the reception 
event and we have ignored atmospheric delays. The actual signals that the GPS satellites transmit 
are continuous-wave circularly polarized radio-frequency signals on two carrier frequencies in the 
L-band centered about: Li ~ 1575.42 MHz and L2 ~ 1227.6 MHz. The carrier frequency signals 
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are modulated by a pseudorandom noise (PRN) code. A GPS receiver makes a phase difference 
measurement, called a pseudorange measurement, which is the phase difference between the PRN 
code received from the satellite and an identical copy of the PRN code that is replicated inside 
the GPS receiver, see Ref.ja] for details of the GPS pseudorange measurement process. The 
pseudorange measurement is essentially made by performing a correlation of the code bits in the 
PRN code received from the satellite with an identical copy of the same PRN code replicated inside 
the GPS receiver. 

Recently, there have been several proposals for synchronizing clocks by making use of entangled 
quantum systems [a, [D, M, M, llil UjJi UM ■ The question naturally arises whether entangled quantum 



systems can be exploited to determine all four space-time coordinates of a user, rather than just 
time. 

In this article, I describe a quantum positioning system (QPS) that is in-principle capable of 
providing a user with all four of his space-time coordinates. This QPS is the quantum analog of 
the classical GPS described above. The QPS is based on entangled photon pairs (biphotons) and 
second order correlations within each pair. The two-photon coincidence counting rate is the basic 
measured quantity. In order to determine his four space-time coordinates, a user of the QPS must 
carry a corner cube reflector, a good clock, and have a two-way classical channel for communication 
with the origin of the reference frame. 

II. INTERFEROMETRIC QUANTUM POSITIONING SYSTEM 

For simplicity of discussion, I assume that space-time is flat with Minkowski 5] coordinates 
(t, X, y, z) and that the user of the QPS is stationary. The complete QPS consists of four biphoton 
sources (entangled photon pairs), four beam splitters and four two-photon coincidence counting 
Hong-Ou-Mandel (HOM) interferometers, see Figure 1. Three of the interferometers are used 
together to solve for the user's position and one interferometer is used to solve for the user's time, 
in a particular reference frame. 

Six spatial points, Rj, where Rj = (xj,yi,Zj), for i = 1,2,3, ..., 6, define the spatial part of 
the reference frame at constant coordinate time t. The six points Rj define three independent 
baselines in pairs, (Ri, R2), (R3, R4), and (R5, Re). The points Rj are assumed to be accurately 
surveyed, so their coordinates are precisely known. Determination of a user's spatial coordinates, 
To = {xo,yo,Zo), is done with respect to this reference frame. A stationary clock in this reference 
frame, say at the origin of coordinates, {x,y,z) = (0,0,0), provides a measure of coordinate time 
t in this 4-dimensional system of coordinates. We neglect all gravitational effects_4] so that the 
user's clock, which keeps proper time r, runs at the same rate as coordinate time t in the reference 
frame defined by the spatial points Rj, so that dr/dt = 1 . Synchronization of the user's clock 
to coordinate time[lj] means that we have a method to compute the integration constant Tq in 
T — To = t. In four dimensional flat space-time, the world lines of the spatial points Rj define a 
tube. Events that are simultaneous in this system of 4-dimensional coordinates occur at constant 
values of coordinate time t, which is a hyperplane that cuts this tube. 

There are several possible modes for the QPS. First, consider a user that must determine all four 
of his space-time coordinates {to, Xo, Vo, Zo)- In a previous paper[l3|, an algorithm has been given 
to synchronize the proper time r on a user's clock to coordinate time, t, without prior knowledge 
of the geometric range from the reference clock (which keeps coordinate time t) to the user clock. 
We assume that this algorithm is employed here to synchronize the user's clock with coordinate 
time t. This algorithm requires a two-way classical channel of communication between the user 
and the reference frame origin, where coordinate time t is kept. The three spatial coordinates are 
determined separately as follows, refer to Figure 1. 



Each baseline, such as the one connected by points, (Ri,R2), contains an entangled photon 
pair (biphoton) source J14l. Ilfj . Ilfil . llTL |l^ located in the baseline. For convenience, we take the 
biphoton source to be at the midpoint of the baseline at point Ei at position ri. Additionally, each 
baseline contains a 50:50 beam splitter and two photon detectors, see Figure |2j For simplicity of 
discussion, we assume that the biphoton source is essentially collocated with the beam splitter and 
two photon detectors at point Ei, see Figure 2. Along the baseline there is a controllable optical 
delay at point Di. The other two baselines contain the same equipment, as shown in Figure 2. 
The QPS works as follows. 

Photon pairs are created at Ei, are sent to positions Ri and R2, and are re-directed to the 
user at the unknown position Tq. The two photon paths are similar, except that one path has 
a controllable optical delay Di. The optical delay is assumed to be calibrated so that we can 
accurately impose an arbitrary delay time[l9|. Next, the entangled photons reflect from the user's 
corner cube reflector at Tq, and return back through the same paths, through points Ri and R2, 
and arrive at a HOM interferometer that is collocated at Ei at position ri , see Figure |2 For 
convenience, we assume that the interferometer is collocated with the photon generation point ri. 
Again, both photon return paths are similar, but one path has the optical delay Di. We have the 
following effective round-trip times for each photon path 

tL = -[|ro-Ri| + |Ri-ri|] (2) 

tR = -[K - R2I + |ri - R2I + (n - 1) d] 

where d is the geometric thickness of the optical delay Di perpendicular to the optical path and n 
is the effective index of refraction for the optical delay Di. The optical delay Di is now adjusted 
until a minimum is observed in the two-photon counting rate at £^i[l^. When the minimum in 
photon coincidence counting rate is observed at interferometer Ei, the effective travel times for 
each photon path are equal, t^ = tR. The interferometer is balanced when the condition ti = Ir is 
satisfied, and a unique minimum is observed in the two-photon counting rate Re- The accuracy with 
which this minimum can be observed depends on the bandwidth Aa; of the band-pass interference 
filters used in front of the photon detectors. 

We get an equation that relates the geometric path lengths to the measured optical delay time 
Ati = (n - l)d/c: 

\ro - Ril -h |Ri - ril = \ro - R2I + |ri - R2I + cAti (3) 

An analogous measurement process is carried out for the other two baselines. For simplicity, I 
assume that the points Ei, E2, and E3 are located at midpoints of their baselines. We then obtain 
the three equations 

|ro-Ri| = |ro-R2| -hsi (4) 

Iro-RsI = |ro -R4I -Fs2 (5) 

Ifo-RsI = Iro-Rel -Fs3 (6) 

where Sj = cAtj for i = 1,2, 3. 



The two-photon coincidence counting rate is given bv|20l.l2ll. |2 



R, = aia2|r?y|2|G(0)|2[l - e'^^'^^^^)'] (7) 

where \V\'^ is the pump intensity in photons per second, ai and 02 are the quantum efficiencies 
of detectors Z?i and D2, r/ is a dimensionless constant and G{t) is the Fourier transform of the 
spectral function (/>, which is the auoto correlation function of the down-converted light 

00 
G{t) = /"(/-(^c^o + uj, ^^0 - u;)e-'^'du; (8) 



The Eqs.Q-® can be solved for the three unknown user spatial coordinates, Fq = (xq. Hoi Zq), in 
terms of the three measured time delays, Ati, At2, Ats, which balanced the three interferometers. 
The measured data consists of photon coincidence count rate vs. optical time delay lengths s,, for 
i = 1,2, 3. Clearly a search must be done of the data to locate the minimum that corresponds to 
equal time of travel along the interferometer arms. The computations can be done at points Ei, E2, 
and E^ . This search to locate the minimum is the quantum analog of the correlation of the PRN 
code in a classical GPS receiver, which was described above. When the three interferometers 
at El, E2, and E^ have been balanced simultaneously, a classical message is sent to the user 
giving him the values of his coordinates r^ = {xo,yo,Zo). Clearly, classical communication is 
needed between the points Ri, R2, and R3 to establish that the interferometers are balanced at 
a given coordinate time t. We imagine that when each interferometer is balanced, a message is 
sent to the origin of coordinates. When three messages are simultaneously received at the origin 
of coordinates (saying that the three interferometers are balanced), Eqs.©-® are solved for the 
user's coordinates Tq = {xq, Ho, Zq) and the user's coordinates are sent to the user through a classical 
channel of communication. 

In the QPS that we describe, there is an apparent asymmetry in the determination of a user's 
spatial coordinates, Tq = {xq, yo, Zq), and in the determination of the user's time. In my view, this 
asymmetry is a reflection of the asymmetric way that space and time enter in the theory of the 
quantized electromagnetic field to give rise to photons as quanta of the field. As mentioned above, 
the time synchronization of the user's clock to coordinate time is done by a method previously 
described by Bahder and Golding_12j. Therefore in what follows, I discuss only the spatial part 
of the QPS. 

With some modification of the above scheme, we may imagine that we could design a similar 
system based on first-order coherence for position determination |20|, |2l|. A single beam from a 
continuous-wave laser can be split and the beams sent on two different paths. However, in such a 
case, there would be an ambiguity that is associated with the wavelength of the light (interference 
fringes will be seen) that is unresolvable in principle. In contrast, in the quantum case (which 
relies on second-order coherence) the ambiguity is resolved because equal propagation times for two 
paths lead to quantum interference: equal travel times for two paths create a unique observable 
minimum in the two-photon coincidence counting rate Re- 

The measured quantities in the QPS are the optical path delays Sj. For a given measured 
value of optical delay, say si, Eq.Q specifies that the user's coordinates must lie on a hyperboloid 
surface with foci at Ri and R2, i.e., a hyperbola of revolution that is symmetric about the baseline 
defined by Ri and R2. The user's position, Tq, is then given by the intersection of three hyperbolas 
given by Eq.Q)-®. Each Eq.Q)-© is just the equation for a baseline in a classical time of arrival 
(TOA) system that records arrival times of classical light pulses (or distinct intensity edges) at 
two spatial reception points Rj. In the case of a classical TOA system, pulse arrival time at four 
locations is needed to determine all four space-time coordinates. In that case, four time difference 



of arrival (TDOA) equations can be formed from four points, and each point is used multiply 
to (effectively) form the baselines. (Taking TDOAs results in a system of three equations where 
the emission event time has cancelled out. ) In the quantum case, since correlations between 
photon pairs are used, we must use three baselines defined by six points Rj, plus an additional 
interferometer for the determination of the user's time. As we will see below, the QPS is an 
interferometric system. 

More fundamentally, and more significant for applications, is that in the classical case of a TOA 
system, we must have good clocks that are synchronized to coordinate time so that accurate pulse 
arrival times at the four Rj reception points can be recorded. A good clock that keeps coordinate 
time is often a difficult requirement to meet in practice 2a]. In contrast, in the quantum case 
two-photon coincidence counts at detectors Diand D2 are measured and only a good "flywheel" 
clock is needed (i.e., a clock having a good short-term stability) to measure photon coincidence 
count rates while the optical time delay is adjusted, to locate the minimum in R^. 

III. GEOMETRIC DILUTION OF PRECISION 

In the case of the classical GPS, the geometrical positions of the GPS satellites determine the 
accuracy of the user's position. This effect is sometimes called the geometric dilution of precision 
(GDOP). We compute the positioning accuracy and the effect of GDOP for the QPS from Eq.@- 
®. These equations give an implicit relation Tq = ro(Ri,R2,R3,R4,R5,R6, si, S2, S3) for the 
user position Tq as a function of the three measured path delays, Sj, and the six baseline endpoints 
Rj. If we knew the error in the measured path length delays, dsi, ds2,and dss, we could compute 
the error in the three components of the user's position vector, dvo = {dXo,dyo,dzo,), from 

1=1 

for constant Rj. However, these errors are statistical in nature, so instead I compute the standard 
deviations ax, (Jy, and o"^, of the user coordinates Xq, Vo, and Zo, as a function of the standard 
deviations o"s^, as2, and ag^, of the measured optical time delays si, S2, and S3. For constant R^ 
for k = 1, ...,6, these standard deviations are related by|2J] 



y \dsij ^1 \ds2J '^ \ds 
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where the partial derivatives are done at constant R^. The lengthy calculation to compute the 
partial derivatives in Ea.()lU|) is done analytically using Mathematica. For simplicity, I assume 
that the error distributions of the Sj are Gaussian and that the three standard deviations are 
equal, a si = (Ts2 = fs3 = cTs- For a spherically symmetric probability distribution of 3-dimensional 
positions r^ = (xq, yo, Zq), the spherical error probable (SEP), which is the radius R within which 
50% of the points lie, is related |25l] to the standard deviations ax = (Ty = (Tz = a hy R= 1.538(T. 
In our case, the probability distribution of Tq is not necessarily spherical. To approximate the 
SEP error metric, we compute a weighted approximation to the SEP metric by defining Rxyz — 
1.538-7=(cr^ + a"^ + a'^)^''^. When the error distribution for Tq is spherically symmetrical, the 



error metrics are equal: Rxyz = R- I consider the effect of GDOP for two cases, one in which 
the interferometer baseUnes are near each other, and the other case where the basehnes are well 
separated, which is the case with classical GPS or a classical TOA system. 

A. Geodetic Positioning System 

First, consider a case where the three baselines coincide with the three Cartesian coordinate 
axes of a reference frame, see Figure Ol Such a case might occur when the baselines are on the 
Earth, and we want to determine the position of an object with respect to a topocentric coordinate 
system. For example, consider the center of the QPS at the origin of Cartesian coordinates and 
an object with a corner reflector at a range of 100m from the QPS, with coordinates (xq. Ho, Zq) = 
(100m)(l, 1, l)/-v/3. Figure 0] shows a plot of contours of constant values of l/Rxyz in the Xo — Vo 
plane at Zq = 100/ v3m, for the interferometer arm (half) length a = 2m and error (standard 
deviation) in optical path ag = 1.0 x 10^^ m. In the contour plot, the position error is Rxyz = 
8.3cm for {xo,yo,Zo) = (100m)(l, 1, l)/\/3, whereas for {xo,yo,Zo) = (30m, 30m, 100/-\/3) the 
error Rxyz = 3.9 cm, which corresponds to the upper right high-accuracy (light-shaded) region in 
FigureEJ On the z-axis at Xq = Ho = and Zq = 100 m the error Rxyz is essentially infinite. Figure 
ISlshows a plot of the error metric Rxyz vs. Xq, for yo = 30 m and Zo = 100m/v3, which corresponds 
to a line in Figure 4 with relatively small error Rxyz- In the high-accuracy light-shaded region 
of Figure EJ for jJq = y^ = 30 m and Zq = 100m/\/3, the dependence of the error Rxyz on the 
baseline length 2a is plotted in Figure also using as = 1.0 x 10~^ m. For a four-meter baseline, 
2a = 4 m, the error is just under 5 cm. 

Note that the position error Rxyz depends linearly on ag, which is the standard deviation (error) 
in measurement of the optical path delay needed to obtain the minimum in two-photon coincidence 
counts Re- The width of this minimum depends on the interference filters in front of the photon 
coincidence counting detectors as well as the pump laser bandwidth [l5|, |2^. Depending on the 
experimental design, this minimum may be measured to better than ag = 1-0 x 10^® m, which was 
used in these plots, and hence accuracies may be better than plotted. 

Finally, I note that the error function Rxyz has a very complex dependence on user coordinates 
{xo,yo,Zo), and as stated earlier, the error function Rxyz also depends critically on the way the 
baselines are distributed, i.e., it depends on the six points R^ for k = 1,...,6, which define the 
baseline endpoints. In the next example, I consider a situation where the baselines do not intersect, 
and thereby the error Rxyz is considerably smaller than for the case considered above, even though 
the distances are larger. 

B. Satellite-based QPS 

Now assume that each point of a baseline, R,, is associated with a different satellite, and that 
the spatial interferometer legs are formed from pairs of points (Ri,R2), (R3,R4), and (R5,R6), 
see FigureEl Specifically, I assume that the points R^, are on low-Earth orbit (LEO) satellites. It 
may seem optimistic that a QPS is feasible with such large baselines because single photons must 
be propagated over these baselines and then reliably detected. However, recently single photons 
have been propagated through the atmosphere and detected over 10 km distance in daylight ,22] , 
and another study concludes that there are no obstacles to create a single-photon quantum key 
distribution system between ground and low-Earth orbiting satellites |28t| . Therefore, a LEO- 
satellite QPS may be possible. 

As an example of the positioning errors in a QPS made from LEO satellites, I take the baseline 
endpoints to be: Ri = (a, -6/2,0), R2 = (a, 6/2,0), R3 = (6/2, a, 0), R4 = (-6/2, a, 0), R5 = 



(-6/(2^/2), -6/(2^/2), a), and Re = (6/(2\/2), 6/(2^/2), a). A plot of this configuration is shown 
in Figure [71 A contour plot of the reciprocal error function, 1/Rxyz, is shown in the Xq — Vo plane 
for Zo = Re/y/S, where Re = 6378 km is the Earth's radius, see Figure |H1 As an example, in the 
calculations below I take the semi-major axis of the LEO satellites to be a = 7360 km and the 
baseline between pairs of satellites as 6 = 20 km. The standard deviation (error) in the measured 
optical delay is taken to be CTs = 1.0 /^m. For a user on the surface of the Earth with coordinates 
{xo,yo,Zo) = (1, 1, 1)-Re/V3 the error is Rxyz = 0.10 cm. For these same parameters, Figure IHl 
shows a plot of the position error Rxyz vs. Xo for Vo = Zq = Re/V^- Note that over a large range 
of Xo-values the error remains below 1cm. Finally, Figure ITUl shows a plot of the position error 
in the radial direction: Rxyz {ro/V^,ro/V^,ro/V^) vs. Tq for the same parameters. On a radial 
line in the (1,1,1) direction, the error remains below 1cm up to Tq = 11680 km. However, near 
1300 km the error rises steeply. This is an example of the complex dependence of Rxyz on user 
position, which was mentioned earlier. 

Clearly, the geometric positioning and layout of the baselines significantly affects the accuracy 
of a user's position. Note that the terrestrial QPS (discussed in the previous section) had a ratio 
of baseline length to user position a/ro = 0.02, whereas this LEO satellite QPS has b/a = 0.003. 
By comparing the baseline layout for the terrestrial QPS and this LEO satellite QPS, it is clear 
that the positioning accuracy is sensitive to the separation and layout of the baselines, but not so 
sensitive to the baseline lengths. Other calculations (not shown) support this conclusion. 

The above calculations for a satellite-based QPS are only meant as an example to illustrate 
the magnitude of errors in position that may be achievable. A significant amount of engineering 
calculations must be performed to design a realist satellite-based QPS. Furthermore, real satellites 
are moving and engineering similar to that used in the classical GPS would have to be done, e.g., 
using Kalman filtering techniques. Obviously, bright sources of entangled photons (biphotons) 
are needed. The calculations above suggest that if properly engineered, a satellite-based QPS may 
achieve position accuracy of objects near the Earth's surface below 1cm. In these calculations, I 
have ignored the time delays introduced by the atmosphere. However, corrections can be made 
for atmospheric effects using multiple colors of photons similarly to what is done with the GPS. 
Perhaps one advantage of the quantum system as compared to the classical GPS is that entangled 
photons exhibit group velocity dispersion cancellation, which may be an important factor for future 
engineering and design of a QPS |29l. I3nl.l3ll.l32||. 

IV. ALTERNATIVE SCENARIOS 
A. Position-Only Determination 

A QPS can be designed to work in several modes, depending on the needs of the user and the 
required scenario. In the above discussion, we have described the case where a user of the QPS 
wants to determine both his spatial and time coordinates, {to,Xo,yo, Zo)- A second alternative is 
that a user may only need to obtain their spatial coordinates, and he may not need the correct 
time. In this latter case, the time synchronization portion of the system is not needed, and the 
user may find their position coordinates {xo, Ho, Zq) by having only a corner reflector and a one-way 
(reception only) classical channel of communication with the reference frame origin, where the 
simultaneity of the three two-photon coincidence counting rate minima is established. 

Another mode of operation of a QPS is where we want to determine the position of an object 
with a corner cube reflector, such as a geostationary satellite. In such a case, information on the 
position of the satellite, Tq = {xo,yo,Zo), is only needed on the ground. The satellite's position 
coordinates can be determined on the ground using a QPS, and only a corner cube reflector is 



needed on the satellite, but no communication channel to the satellite and no on-board clock is 
needed. 

B. User Carries QPS Receiver 

The scenarios that we have described above are ones where the measurements (adjusting the 
optical delays) and the calculations (to compute Tq) are done near the origin of the reference 
frame. In a classical GPS receiver, the computations (correlations of PRN codes to at least four 
satellites) are done locally in the user's GPS receiver that the user carries with him. The QPS 
analog of this classical GPS scenario is a setup where the biphotons are generated at points Ei^ 
E2, and E^, but the user carries with him the 50:50 beam splitters and photon detectors. In this 
scenario, the user controls (and carries with him) the optical delays, see Figure ^ and he locally 
measures the optical delays si^ 52, and S3. The user must receive a classical message consisting 
of the coordinates of baseline endpoints. Hi, i = 1,...,6, and then he must solve the Eq.©-® 
for his position Tq = {xo,yo,Zo)- In such a case, there are no clocks on-board the broadcasting 
satellites (located at positions Rj), however, the user must carry a clock with short term stability 
to determine rate of photon coincidence counts from each of the three baselines (associated with 
spatial positioning) and also he must do coincidence counting for time determination (if time is 
needed). For the three spatial baselines, optical propagation is then one-way (using the satellite 
positions as a primary reference system, see below) from satellites to QPS user receiver. For 
time synchronization, however, as mentioned previously, the optical propagation must be two-way 
(when using the method of Bahder and Golding). In essence, for each of the four channels, the QPS 
receiver consists of a beam splitter, two single-photon detectors, and a controllable optical delay. 
All four space-time coordinates can be obtained by a user in this way. One clock in the reference 
frame must have long-term stability to define coordinate time, and another clock in the QPS user 
receiver can have short-term stability. Note that the satellites do not need to carry clocks, because 
their positions can be used to define the primary system of coordinates. This type of QPS is a 
close analog of the classical GPS. 

V. QPS SPACE-TIME COORDINATES 

The satellites at baseline points Rj can be taken to define the primary system of reference, 
even though the points Rj change with time. The quantities measured by a user of such a QPS 
are then {sq, 81,32,83), where sq is the optical time delay (in the HOM interferometer) that will 
provide the user with coordinate time in this coordinate system (using the Bahder and Golding 
method), and (si, 52,^3) are the three optical delays in the three interferometers for position 
determination. The quantities (sq, si, S2, S3) are then to be regarded as generalized 4-dimensional 
space-time coordinatespj], sq is a time-like coordinate and {81,82,83) are space-like coordinates. 
Within the context of general relativity, such coordinates are as good as any other coordinates, 
and they enter into the metric c^dr'^ = gijdsidsj of the flat space-time assumed in this work. Of 
course, a transformation from the QPS space-time coordinates, {sq, 81,82, S3), to an Earth-centered 
inertial (ECI) system of coordinates, say (t,x,y,z), is of interest for astrodynamic applications. 
Such a transformation can be done approximately by conventional means of tracking the satellites 
(at baseline points Rj). 

It is interesting to remark that the QPS allows the direct measurement of 4-dimensional space- 
time coordinates. Previously, it was believed that space-time coordinates were not measurable 
quantitiesSaSa. Of course, the QPS coordinates (.0, .1, .2, .3) are real physical measure- 
ments, and it is well-known that real measurements are space-time invariants under generalized 



coordinate transformations 1331. 



VI. SUMMARY 

I have presented a conceptual scheme for an interferometric quantum positioning system (QPS) 
based on second order quantum coherence of entangled photon pairs (biphotons). A user's spatial 
coordinates are determined by locating three unique minima in three different two-photon counting 
rates, associated with three HOM interferometers built on independent baselines. The spatial 
portion of the QPS is similar to a classical TOA system, however, a classical TOA system requires 
synchronized clocks that keep coordinate time, which is often a difficult requirement to meet. 
In contrast, the QPS only requires a clock having a short-term stability to measure two-photon 
coincidence counting rates while the optical time delay is adjusted (to locate the minima in the 
two-photon coincidence counting rate Re). Bright sources of entangled photons (biphotons) are 
needed. 

Several different scenarios were considered for a QPS: one is a terrestrial system and the another 
is space-based. In both cases, I computed the accuracy of a user's position as a function user 
position. The function that describes the errors in position has a complex spatial dependence. In 
the case of the terrestrial QPS, the position accuracy was relatively poor because the baselines 
were located near each other. This could be dramatically improved by moving apart the baselines. 

As an example of a satellite-based QPS, I have proposed a LEO-satellite QPS. Neglecting 
atmospheric effects, calculations suggest that the position accuracy Rxyz of such a QPS can be 
below the 1 cm-level for an error (standard deviation) in the optical delays ag = 1.0 ^m associated 
with the minima in two-photon counting rates Re- The complex dependence of Rxyz on user 
position suggests that significant engineering must be done to design a realistic QPS. 
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FIG. 1: One baseline is shown for the quantum positioning system (QPS). Points Ri and R2 on the baseline 
contribute to the definition of the reference frame for spatial positioning. Box Ei contains an entangled 
photon (biphoton) source and 50:50 beam spliter, see Figure 2. The quantity Di is a controllable, calibrated 
optical delay. 
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FIG. 2: An expanded view of the contents of each of the three boxes Ei, E2, and E^, which are located on 
the three basehnes, one of which is shown in Figure 1. Each box contains an entangled photon (biphoton) 
source, a 50:50 beam splitter, and two single-photon detectors -Diand D2, to perforin photon coincidence 
counting. 
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FIG. 3: The baselines are shown for a possible terrestrial QPS that might be used on the Earth, 
baselines lie along the x, y, and z axes, are of length 2a and are orthogonal to each other. 
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FIG. 4: A plot of the contours of constant 1/Rxyz is shown in the Xo — yo plane at Zo — 100/\/3i 
Light-shaded areas are small values of Rxyz- Units on both axes arc meters. 
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FIG. 5: The error Rxyz vs. Xo is plotted for yo = 30m and Zo = lOOm/vo, which corresponds to a hnc in 
Figure 4 with relativily smaU error Rxyz- The same parameters are used in this plot as in Figure 4. Units 
on both axes are meters. 
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FIG. 6: The error Rxyz vs. a (half the interferometer basehne length) is plotted for Xo ^ Vo = 30 m, and 
Zo ~ 100m/-\/3, which corresponds to the high-accuracy light-shaded region in upper right of Figure 4. 
Units on both axes are meters. 



17 




z 



FIG. 7: A schematic of the LEO satellite QPS is shown. Pairs of satellites orbiting Earth, shown by 
connecting lines, form the interferometer baselines of length b. Example numbers used in this calculation 
have baseline 6 = 20 km and LEO satellite semi-major axis a = 7360 km . 
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FIG. 8: Contours of constant reciprocal position error, 1/Rxyz, are shown in the Xo — yo plane for Zq = i?e/V3 
and (7s = 1.0 /xm. Lighter-shaded areas are smaller values of error Rxyz- The semi-major axis of the LEO 
satellites is taken to be a = 7360 km and the baselines (satellite pair separation) 6 = 20 km. 
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FIG. 9: Plot of position error, Rxyz vs. Xq, shown for yo 



Re/v3 with (jg = 1.0 /ini. The semi-major 



axis of the LEO sateUites is taken to be a = 7360 km and the basehnes (sateUite pair separation) b = 20 km. 
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FIG. 10: Plot of position error in the radial direction, Rxyz vs. 7 

as = 1.0 /im. The semi-major axis of the LEO satellites is taken to be a 

(satellite pair separation) 6 = 20 km. 



where Tq — ya^^+l/J+^J and 
7360 km and the baselines 



